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We show that a simplified two-band model describing toroidal magnetic order in two-dimensional 
crystal is entirely equivalent to the well-known Haldane model of a honeycomb lattice in periodic 
internal magnetic field with zero total flux through the unit cell. Crystals with toroidal ordering can 
be considered as real physical systems which should exhibit topologically nontrivial state, predicted 
by Haldane, that is similar to the quantum Hall effect regime without any external magnetic field. 
Analogous description is formulated for quantum spin Hall effect without any spin-orbit interaction. 
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PACS numbers: 73.43.-f, 75.85.+t, 74.90.+n 

Highly speculative model of a two-dimensional (2D) 
honeycomb lattice with two nonequivalent sub-lattices in 
periodic internal magnetic field, introduced by Haldane 
more than twenty years ago [ljj, directly predicts a new 
state of matter similar to the integer-valued Quantum 
Hall Effect (QHE) regime without any external magnetic 
field. It is presupposed in this Haldane model that in- 
ternal magnetic-flux density with full symmetry of 2D 
lattice results in zero total flux through the unit cell. In 
the framework of the Haldane model, time-reversal in- 
variance and 2D space-inversion symmetry are broken by 
the internal magnetic field and the non-equivalence of the 
sublattices, respectively. Tight-binding model [lj takes 
into account hopping matrix elements between nearest 
and second-neighbor sites (the nearest sites of a sublat- 
tice) and the periodic vector potential A(r) of the in- 
ternal field is included into this model through the phase 
multipliers exp [i(e/hc) J Adr] of the hopping matrix ele- 
ments. Thus, magnetic field in the model [l[ becomes ap- 
parent through the vector potential in the phase factors 
of the wave functions similar to that in Bohm-Aharonov 
effect. Since internal magnetic field breaks time-reversal 
invariance, quasiparticle spectrum of the system with 
broken spatial inversion turns out to be nonsymmctric 
with respect to inversion of quasi- momentum, k — k. 
In addition, edge chiral fermion excitations, typical of the 
QHE regime, arise without partners of opposite chirality. 

It should be noted that, simultaneously with and in- 
dependently of the paper Ref . [l[ , Volovik @, H| has pre- 
dicted and studied in detail the half-integer QHE in the 
A-phasc of liquid helium-3 without both external mag- 
netic field and rotation of the container. 

In the QHE regime, quantizing external magnetic field 
normal to the plane of 2D crystal both transforms 2D 
metal into 2D insulator in which degenerated Landau 
levels play role of the bands of allowed states and also 
determine a direction of a drift of the cyclotron orbit 
center along the boundary of the crystal. Finite distance 
between the neighboring Landau levels and the absence 
of the backscattering in the edge chiral states result in 
a stability of the ground state of the QHE regime with 
respect to relatively weak structure perturbations. Com- 



paratively weak internal periodic magnetic field in the 
Haldane model breaks time-reversal symmetry but in it- 
self cannot create an insulating state in the form of a 
system of Landau levels like that in the QHE regime. 
However, as was also proposed by Haldane some- 
times the internal field could be associated with a mag- 
netic ordered state with broken time-reversal invariance 
that might result in an insulating gap in the quasiparticle 
spectrum. Finite insulating gap and absence of backscat- 
tering could lead to a robust Haldane state similar to the 
QHE regime. 

Thouless et al. [H have shown that QHE regime dif- 
fers from the conventional band insulating state by the 
integer-valued topological invariant: in the case of con- 
ventional (trivial) insulator, such a topological invariant 
turns out to be zero, n — and edge states are absent; 
ri 7^ may be referred to a topological insulator. 

Edge states can also arise in the topological insulator 
with spin-orbit interaction which does not break time- 
reversal symmetry 5, 6] . Such edge states correspond to 
a transfer not of electric charge but of spin. Electrons 
with opposite spin projections move in opposite direc- 
tions along crystal boundaries. This becomes apparent 
as spin quantum Hall effect [!, 0] . 

As it was pointed out by Haldane [l[ and Volovik 
3] , topologically nontrivial insulating state with electric 
charge transfer over the edge states could be realized in 
2D systems with magnetic order breaking time-reversal 
invariance of the Hamiltonian. However, one can con- 
clude that such states have never been experimentally 
realized till now. In this Letter, we show that two-band 
model of excitonic insulator Q with electron-hole pair- 
ing leading to toroidal magnetic order turns out to be 
equivalent to the Haldane model after very natural sim- 
plification. We believe the Haldane state can be realized 
just in 2D insulating crystals with toroidal magnetic or- 
der. 

Such insulating crystals with spontaneous toroidal mo- 
ment of the unit cell in the ordered phase are investigated 
theoretically long ago [1]. Since then, a great number of 
crystals and artificial 2D heterostructures with toroidal 
orderring were discovered Q mainly because of their re- 
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lationship to the magnetoelectric effect. 

Toroidal magnetic moment t arises in the multipole 
expansion of the electrodynamic vector potential A(r) 
beginning with the second order term which can be pre- 
sented as a sum of the contributions from the magnetic 
quadrupole and toroidal moments [l(| • The toroidal part 
of the vector potential created by orbital currents at some 
distant point r can be written as Q 

A(r) = V(tV)r" 1 + 47r-r(r), (1) 

where 

r(r) = — [r X [r X j(r)}} (2) 

is the toroidal moment density, j(r ) is the orbital electric 
current density. Nonzero toroidal moment corresponds to 
a special distribution of orbital currents in Eq. ([2]). As a 
classical example, one can consider the so-called poloidal 
currents flowing through a solenoid that is bent into a 
torus. 

One can imagine that all orbital currents are shrunk 
into the origin of coordinates [9(. Then the toroidal mo- 
ment density ([T]) should be written as t — > t ■ 5(r). The 
first term in Eq. ([TJ can be eliminated by appropriate 
gauge transformation so that the vector potential takes 
the form A(r) = 47rr(r). 

It should be noted that toroidal moment can be created 
not only by persistent orbital currents but also by an 
appropriate configuration of magnetic dipole moments, 
such as head-to-tail arrangement of spins [9( . 

One can see that the toroidal contributions into both 
vector potential and magnetic flux density B = rotA 
turn out to be localized inside a bounded domain of the 
real space in which toroidal moment density is nonzero. 
In the case of crystalline solid, such a domain is the unit 
cell so that both vector potential and magnetic flux den- 
sity should be periodic functions with basic periods of 
the solid similar to that proposed in the Haldane model. 

Since toroidal moment density t(t) is £-odd polar (r- 
odd) vector, the mean field order parameter originat- 
ing from toroidal ordering should change sign upon both 
time reversal and spatial inversion. As direct product 
of electric and magnetic field strengthes has the same 
symmetry toroidal order may be apparent as a linear 
magneto-electric effect. Corresponding contribution into 
the free energy of the form F me = —a ik E i H k where 
a ik is a magneto-electric tensor. Magnetic field induced 
electric polarization and electric field induced magneti- 
zation are determined by this tensor as P i = cn ik H k and 
M i = a ki E k , respectively. 

From the symmetry consideration, magneto-electricity 
should be one of the distinctive properties of the toroidal 
ordered state fy: P = — [txH], M = [txE]. Magneto- 
electric tensor can be decomposed into a pseudoscalar, 
a polar vector and a symmetric traceless tensor. Thus, 



non-diagonal components of the magneto-electric tensor 
should be nonzero and antisymmetric, a ik = — a ki [111 ]. 

The bulk adiabatic magneto-electric tensor determines 
the surface Hall impedance in a magneto-electric crystal 
cn^^n^ where c is light velocity. This 
thermodynamic relation is used in Ref. [l[ to calculate 
quantized Hall conductivity cr xy . 

A rise of the toroidal ordering can be a result of an elec- 
tron phase transition. It is convenient to describe such 
a transition in the framework of the model of the exci- 
tonic insulator with high-temperature (more symmetric) 
phase being 2D semimetal with equally centered isotropic 
electron and hole bands [1]. 

Bloch Hamiltonian of such two-band semimetal can be 
written as 

H p = e + (p) r°+e_ (p) f 3 + ( 7 p) f 2 , (3) 

where f° is the unit matrix, f 1,2 ' 3 are the Pauli matrixes, 
2e±(p) = e 1 (p)±e 2 (p), e s (p) is electron dispersion in the 
band s (electron energy is counted from the chemical po- 
tential), 7 = Pi 2 /m is a real vector, p 12 is interband ma- 
trix element of momentum, m is electron mass. Note that 
p 12 =/= when, in particular, the bands of the semimetal 
are formed by the hybrid states of different parity, for 
example, s and p states. 

Fermi contour (FC), that is a boundary separating fill- 
ing and vacant states in the electron and hole bands, is 
determined by the Fermi momentum p F . In the case of 
isotropic bands the FC is a circle with radius p F . One 
can choose praxis along the vector -y. For low-energy 
elementary excitations a deviation k of quasimomentum 
from the Fermi momentum p F is small. Therefore, in the 
first approximation, one can replace the quasimomentum 
p with p F in the non-diagonal terms of Hamiltonian (|3]) 
whereas the diagonal terms depend on k directly because 
energies s 1 (p) and e 2 (p) are counted from the chemical 
potential. 

Hamiltonian © is characterized by gapless anisotropic 
quasiparticle spectrum 

i? li2 (fc) = e + (fc)± v /e 2 _(fc) + (7p F ) 2 , (4) 

which is symmetrical with respect to quasimomentum 
inversion, k — > — k. 

To obtain both time-reversal and space-inversion 
symmetries to be broken one can consider the low- 
temperature ordered phase originating from electron-hole 
pairing resulting in nonzero toroidal moment. The or- 
der parameter is determined by an anomalous average 
(^ipa^2pa') corresponding to singlet (er = a') or triplet 
pairing of electron and hole belonging to different bands. 
Here, c spa annihilates electron with quasimomentum p 
and spin a in the band s = 1,2. 

It is well known that orbital and spin current states 
corresponding to singlet and triplet electron-hole pairing, 
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respectively, can result in the ordered states with pure 
imaginal order parameter [l3j]. 

We restrict ourselves to the case of singlet electron- 
hole pairing resulting in a pure imaginal order parameter 
iA(fc) where A(fe) is a real function of quasimomentum 
k. In such a case, toroidal order arises from poloidal 
charge currents. To associate the order parameter with 
the toroidal moment one should use the relation t = 7A 
following from the interrelation between the char ge c ur- 
rent density and toroidal moment of the unit cell 1J] . 

The singularity in the electron-hole pairing channel 
that could lead to the insulating ordered state at low 
temperature can arise under nesting of electron disper- 
sion, e 1 (k) = —e 2 (k) = e(k). We suppose that such a 
condition is fulfilled in a vicinity of the FC. Then the 
mean field Bloch Hamiltonian corresponding to singlet 
electron-hole pairing with a rise of a toroidal magnetic 
moment of the unit cell and pure imaginal order param- 
eter A(fc) takes the form 

H k = h(v F k)T 3 + [( 1PF )+A(k)}f 2 , (5) 

where electron dispersion is presented by linear function 
of quasimomentum, s(k) = Ti(v F k). Here, v F is Fermi 
velocity. Later on, we restrict ourselves to the case of 
isotropic (s-wave) toroidal order parameter A = const. 

The spectrum of the high-temperature phase, when 
A = 0, 

E h2 (k) = ±^h 2 {v F ky + { 1PF )\ (6) 

turns out to be gapless because of the fact that jp F = 
for two antipodal directions of the p 2 -a,xis. Two points on 
the FC corresponding to such directions form two Dirac 
cones with linear dispersion at small k. 

In the low-temperature phase, when A =/= 0, quasipar- 
ticle spectrum corresponding to Bloch Hamiltonian §5§ 
becomes gapped, 

E 1>2 (k) = ±^h 2 {v F k) 2 + { 1 p F + A) 2 . (7) 

Similar to the Haldane model |l| , the spectrum turns out 
to be asymmetric with respect to k — > —k: E 1 2 (— fe) 7^ 
E 1:2 (k). 

Note that quasiparticle dispersion (J7]) depends on k not 
only over the first term under the square root but also 
over the scalar product depending on quasimomentum 
polar angle ip. Since vector 7 is directed along p 1 axis, 
{"IPf) > at \ip\ < 7r/2; in the opposite case, {~fP F ) < 0. 

The energy gap in the spectrum ([7]) depends on quasi- 
momentum considerably. In particular, it equals to 2A 
at the points of intersection of the FC and p 2 -axis. In the 
case of ~YP F > 0, the gap monotonically increases from 
2A when polar angle ip between 7 and p F varies from 
±7r/2 (these points corresponds to jp F — 0) to zero (if 
tp = 0, the gap equals to 2y / "/ 2 p F + A). 



If ('YPp) < 0, the energy gap is a nonmonotone func- 
tion of the polar angle. It vanishes at the points of the FC 
with polar angle <p — ±(fi c corresponding to "fP F + A = 0. 
Thus, if A is small enough the spectrum (J7J remains gap- 
less. The gap can arise in those directions of the momen- 
tum space where the toroidal order parameter exceeds 

1P F '- A > IPf- 

In the case A < jp F , one can expect that Bloch Hamil- 
tonian ([5]) should be related to a semimetal with in- 
sulating gap on some part of the FC corresponding to 
— 1 < cosi^ < —A/"/p F . Distribution of semimetal and 
insulating constituents in the real space corresponding 
to a peculiar intermediate state may be highly compli- 
cated. One can think that the phase transition between 
semimetal and insulating phases should be occurred un- 
der the condition that the insulating toroidal order pa- 
rameter exceeds several critical value, A > A c , where 
A c can be considered as a minimum value of the toroidal 
order parameter. 

Thus, owing to ^-dependence, the mean field Bloch 
Hamiltonian §5§ describes rather complicated low- 
temperature insulating phase with toroidal magnetic or- 
der. Hamiltonian ([5]) can be simplified considerably if 
one replaces scalar product ~/p F with constants corre- 
sponding to the averages over the polar angle in the lim- 
its of the domains of constant sign of this scalar product, 
jp F -> aj c p F . Here a = sgn(jp F ), j c = A c /p F . 

In just the same way as in the case of the Haldane 
model, mean field Bloch Hamiltonian <j5j) can be reduced 
to two independent effective Hamiltonians, 

H ak = h(v F k)f 3 + [a lc p F + A]f 2 , (8) 

which are distinguished by discrete parameter a = ±1, 
which connects the terms breaking space inversion and 
time reversal symmetries ("YPp) + A — > a^ c p F + A. 

One can see a direct analogy between the pairs of ef- 
fective Bloch Hamiltonians corresponding to extremely 
simplified problem of the excitonic insulator with toroidal 
magnetic order |8]) and the Haldane model |l| described 
by the effective Hamiltonians (we use somewhat other 
designations with respect to the original paper [H) 

#iS? = Kv F (k x T 2 - kyT 1 ) + m a v 2 F T 3 (9) 

which can be considered as a presentation of a traceless 
2x2 matrices in the form of a linear combination of 
the Pauli matrices. The two effective Hamiltonians (jHJ 
coincide with ([9|) nearly verbally if one considers more 
general case of the ordering with nonzero both real and 
imaginal parts of the order parameter, A = A' + iA" = 
\A \ e 1 ^. Here, <j> is a phase of the order parameter (4> — 
7r/2 in the case of pure imaginal A). 

The antipodal directions of quasimomentum, differing 
with sign of a, turn out to be considerably nonequivalent 
at A ^ 0. Energy gap E g+ — 2(A C + A) is an increasing 
function of A in the case of a direction corresponding to 
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a = +1, whereas the gap E g _ = 2| A c — A| for a direction 
corresponding to a = — 1 exhibits a minimum at A = A c 
in which E g _ = 0. Change of sign of the gap corresponds 
to the band inversion similar to that in the systems with 
quantum spin Hall effect It is important that band 
inversion in excitonic insulator is not a result of the spin- 
orbit interaction so that one can expect that the Haldane 
state can be observed at rather high temperature. 

Approximate Hamiltonians H ak corresponding to two- 
band model of excitonic insulator with toroidal magnetic 
order and effective Hamiltonians of the Haldane 

model describing 2D honeycomb lattice are not identi- 
cal. Nevertheless, both and H ak result in the same 
structure of quasiparticle spectrum. In the case of exci- 
tonic insulator we have 



Unitary matrix connecting two-band operators (|5J) and 
(|9]) depend on three angle variables, 9, (p. Matrix ele- 
ments can be presented in the form = U£ 



7?9 

e cosy), 



EUk) = ±Jh\v F k)* + (aA c + A)2 



(10) 



where the second term under the square root plays role 
of the mass term m a defined in the Haldane model [l| as 

m a = M — 3v3at 2 sin0 — > aA c + |A|sin^>; 

here ±M are the energies of the sites of two different 
sublattices of the honeycomb lattice, t 2 is a matrix ele- 
ment corresponding to the hopping between the second 
neighbor sites, <j> is an additional phase corresponding to 
such a hopping in the Haldane model [l[. 

Note that triplet electron-hole pairing described by 
pure imaginal order parameter can be associated with 
poloidal spin currents in the volume [15| and therefore 
should be apparent in spin quantum Hall effect. Triplet 
pairing can be described approximately by the Hamilto- 
nian in the form of a direct sum of two pairs of effective 
Hamiltonians H kcra , similar to © or ©, for different 
spin projections a, where a = ±1 for a =*f whereas 
a — =pl for a Similar model was introduced by 

Bernevig, Huges and Zhang [6] to study the quantum 
spin Hall effect in HgTe/CdTe quantum wells. 

Any Hermitian operator is fully determined by its 
spectrum E n and a complete set of eigenfunctions tp nv , 
where quantum number v takes into account a degener- 
acy of the eigenvalue E n . Such a set can be used as an 
orthonormal basis of the Hilbert space in which operator 
H^- ' acts. Let us consider another Hermitian operator 
H acting in the same space and suppose that and 
H have the same spectrum E n . This fact indicates that 
the two operators are equivalent, 



H = UH ( - 0) U\ 



(11) 



where U is a unitary operator. Indeed, a unitary trans- 
formation preserves the spectrum in a new orthonormal 
basis Uip ' . 



U 12 = —U21 = e sin<y9. The explicit form of these matrix 
elements can be directly obtained from Eq. (fTTj) . 

One can see that the excitonic insulator with toroidal 
magnetic order described by the approximate Bloch 
Hamiltonian (JS) and the Haldane model [1| described by 
the Bloch Hamiltonian ([9]) are equivalent that is can be 
related to a class of models in which the Haldane state 
can be apparent. 

To create the Haldane state one needs only the time- 
reversal symmetry break in 2D crystal. Therefore, it 
seems that artificial heterostructures in the form of quan- 
tum well systems may be the materials in which such a 
state can be discovered. To obtain desired properties of 
the electron spectrum of the heterostructure one can eas- 
ily correct this spectrum by external in-plane magnetic 
field [IH or by varying of magnetic impurity concentra- 
tion. We believe that precisely 2D quantum wells with 
toroidal magnetic order are the most suitable systems to 
realize the Haldane state. 
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